Josephson Effect in Fulde-Ferrell-Larkin-Ovchinnikov Superconductors 



&\ 
On 

o 

Q 

o 

(N 



Kun Yang and D. F. Agterberg 
National High Magnetic Field Laboratory and Department of Physics, Florida State University, Tallahassee, Florida 32310 

(February 1, 2008) 

Due to the difference in the momenta of the superconducting order parameters, the Josephson current 
in a Josephson junction between a Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) superconductor and 
a conventional BCS superconductor is suppressed. We show that the Josephson current may be 
recovered by applying a magnetic field in the junction. The field strength and direction at which 
the supercurrent recovery occurs depend upon the momentum and structure of the order parameter 
in the FFLO state. Thus the Josephson effect provides an unambiguous way to detect the existence 
of an FFLO state, and to measure the momentum of the order parameter. 
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It was suggested more than thirty years ago by Fuldc 
and Ferrell [[[), and Larkin and Ovchinnikov Q, that an 
inhomogeneous superconductor with an order parameter 
that oscillates spatially may be stabilized by a large ex- 
ternal magnetic or internal exchange field. Such a Fuldc- 
Ferrell-Larkin-Ovchinnikov (FFLO) state has never been 
observed in conventional low-T c superconductors. Re- 
cently it has attracted renewed interest in the context of 
organic, heavy-fermion, and high-T c cuprate supercon- 
ductors |^jl9| ■ These new classes of superconductors are 
believed to provide conditions that are favorable to the 
formation of FFLO state, because many of them are i) 
strongly type II superconductors so that the upper criti- 
cal field H C 2 can easily approach the Pauli paramagnetic 
limit; and (ii) layered compounds so that when a mag- 
netic field is applied parallel to the conducting plane, the 
orbital effect is minimal, and the Zeeman effect (which is 
the driving force for the formation of FFLO state) domi- 
nates the physics. Indeed, some experimental indications 
of the existence of the FFLO state have been reported 

The main difficulty in the experimental search for the 
FFLO state is that just like the BCS state, the FFLO 
state is a superconducting state. The distinction be- 
tween the two is a subtle difference in the structure of 
the superconducting order parameter, which is difficult 
to detect using ordinary experimental techniques. Pre- 
vious experiments have focused on thermodynamic sig- 
natures of possible phase transitions from the BCS to 
FFLO state, which is believed to be first order (see, how- 
ever, Ref. (|]). But such signatures can be caused by 
other phase transitions that have nothing to do with the 
superconducting order parameter. Thus it is very diffi- 
cult to establish the presence of an FFLO state this way 
without any ambiguity. 

In this paper we propose using the Josephson effect 
to detect the existence of FFLO states. Our proposal 
has some similarity in spirit to the basic ideas behind 
the "phase sensitive" experiments [^0| that established 
the predominant d-wave symmetry of the order parame- 



ter in cuprate superconductors. Specifically, we predict: 
(i) The Josephson current in a Josephson junction be- 
tween a conventional BCS superconductor and an FFLO 
superconductor is suppressed, due to the difference in 
momenta of the order parameters, (ii) The Josephson 
current may be recovered by applying a properly chosen 
magnetic field in the junction, with field strength and 
direction depending on the momentum of the order pa- 
rameter of the FFLO superconductor; it thus provides a 
way to measure the momentum of the order parameter 
directly. 

In the rest of the paper, we demonstrate the above ef- 
fects first by using the Ginsburg-Landau theory, and then 
by presenting a microscopic derivation. We discuss the 
experimental implication and feasibility of our proposal 
toward the end of the paper. 

Ginsburg-Landau Theory. The effects we predict are 
most easily demonstrated using a Ginsburg-Landau the- 
ory. For simplicity we consider a two-dimensional BCS 
superconductor, described by a spatially dependent su- 
perconducting order parameter ^bcs{ v )i which is cou- 
pled to a two-dimensional FFLO superconductor, de- 
scribed by an order parameter ^fflo{ v ) Pi]] - We con- 
sider the two Josephson junction geometries shown in 
Figure 1. Since the physics for the two geometries is sim- 
ilar we focus our discussion on geometry a) and simply 
state the results for geometry b). For geometry a) the 
Josephson coupling term in the free energy takes the form 
(in the absence of any magnetic field) 



Hj = -tJ d 2 v[m FFLO (v)^ BCS 



(r) + c.c. 



(1) 



In the ground state of a BCS superconductor, \P Bcs{ r ) = 
■00 is a constant. However, in an FFLO superconductor 
the order parameter is a superposition of components 
carrying finite momenta: 



(2) 



and is oscillatory in space. In the absence of magnetic 
flux inside the junction, the total Josephson current is 



f 



t J d 2 r** BCS {r)* FFLO (r) 



(3) 



Clearly, due to the oscillatory nature of the integrand, 
the Josephson current is suppressed in such a junction. 
The same result is clearly true for geometry b). 

Mathematically, the reason that the Josephson current 
is suppressed here is similar to the suppression of Joseph- 
son current by an applied magnetic field in an ordinary 
Josephson junction between two BCS superconductors 
p2| . However, the physics is very different: here the 
suppression is due to the spatial oscillation of the order 
'parameter in the FFLO state, while in the case of ordi- 
nary Josephson junction in a magnetic field, the phase of 
the Josephson tunneling matrix element is oscillatory (in 
a proper gauge choice). Nevertheless, the mathematical 
similarity allows these two effects to cancel each other 
and restore the Josephson current, as we demonstrate 
below. 

Consider geometry a). Imagine applying a parallel 
magnetic field B_!_z parallel to the planes, where z is a 
unit vector along the normal direction of the plane. Using 
the gauge A(r) = r x B we have A(r) = A(r)z. The ap- 
pearance of the A field affects the phase of the Josephson 
tunneling matrix element only; the in-plane properties 
in the two individual superconductors are unaffected be- 
cause the A field is perpendicular to the planes. Specif- 
ically, the Josephson coupling term in the free energy 
becomes 



Hj = -t 



„ 2ieA(r)d 

cTr[e ~- * 



FFLO 



(r)$BCs(r)+c.c], (4) 



and the total Josephson current takes the form 



I j = Im 



d 2 



re" 



eA(r)d 



**BCs(r)*FFLo(r) 
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Here d is the distance between the two planes. It is clear 
from the equation above that the two oscillatory factors 
cancel each other when 



B 



he . 

2ed* 



x k r . 



(6) 



At this particular B, the Josephson current gets partially 
restored. One can determine the momenta of the order 
parameter in the FFLO state by searching for the B's 
that restore the Josephson current. 

For geometry b) consider applying a perpendicular 
magnetic field B || z. Arguments similar to those of 
the last paragraph imply that the Josephson current be- 
comes restored only when the momentum of the FFLO 



order parameter is along the junction (note that typically 
this momentum will lie along high symmetry directions of 
the superconductor). In this case the Josephson current 
is partially restored when 



B 



he 



2e(d + Ai + A 2 ) 



(7) 



where Ai (A 2 ) is the penetration depth of the FFLO 
(BCS) superconductor for the field applied normal to 
the superconducting plane. The appearance of the fac- 
tor d + Ai + A 2 in Eq. ^ is because the magnetic field is 
assumed to enter each superconductor a distance equal 
to the penetration depth from the edge of the junction. 

Microscopic Derivation. These effects can also be 
demonstrated using a microscopic approach, as we il- 
lustrate below. For simplicity and clarity, we take 
the zero temperature limit, and initially consider the 
usual Josephson effect between two BCS superconduc- 
tors. When two superconductors are placed in proximity 
so that electrons can tunnel from one to the other, the 
total energy of the system contains a term that depends 
on the phase difference 4> between the order parameters 
of the two superconductors: 



Ej(4>) = E j cos( 
and the Josephson current is 

T 2eEj ■ 

Ij{<P) = — r— sin 



(8) 



(9) 



Thus the calculation of Ij((j>) is equivalent to the cal- 
culation of the Josephson energy Ej(4>), which we per- 
form below. We use the following tunneling Hamiltonian 
MM: 



H 7 



pqcr 



(r P , 



h.c), 



(10) 



where p and q label single electron momentum eigen- 
states in the two different superconductors, and a is the 
spin label (this tunneling Hamiltonian can be used for 
both geometries since only the matrix elements T p q will 
differ). The mean field Hamiltonian for a BCS supercon- 
ductor (s-wave) takes the form: 

per p 

= ^ £ P (QpQp +/3 P /3 P ) + const., (11) 
p 

where A; is the order parameter of the lower super- 
conductor, a and /3 are Bogliubov quasiparticle opera- 
tors, and E p = ^Je p + |A;| 2 is the quasiparticle energy. 

When both lower and upper superconductors are BCS, 
a second-order perturbation calculation in Ht gives rise 
to the following expression for the Josephson energy (we 
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have neglected terms in the second-order perturbation 
that are independent of the phase difference (f>): 



T A* A 



EpEq(E p + Eq 



+ c.c. 



(12) 



In general the tunneling matrix elements T p q are com- 
plex. A fact that is crucial to the existence of the Joseph- 
son effect, as originally emphasized by Josephson himself 
|p3f , is that in the presence of time-reversal symmetry, 
T_ p ._ q = T* ; thus there is no phase oscillation in the 
summation of Eq. (12). Breaking time-reversal symme- 
try (e.g. by applying an external magnetic field) tends to 
suppress the Josephson effect by introducing oscillatory 
phases to the terms in Eq. (12). 

We now turn the discussion to the case where the upper 
superconductor is in an FFLO state, while the lower one 
is still a BCS superconductor. We first consider zero ex- 
ternal magnetic field, and for simplicity assume the order 
parameter of the FFLO state has a single momentum k. 
The mean-field Hamiltonian in this case takes the form 

§§: 

H u = £ qC qCT c q(T + E z ^( c qT c iT - c qi c ql) 

qcr q 
+ X/(^ uC q+k/2T C -q+k/21 + ^' C -) 

q 

= $}(-Eq + z qK"q + ( E i - ZJPtPi] + const -> ( 13 ) 
q 

where Ez is the Zeeman splitting between the up and 
down spin electrons, and Zq = Ez + vpQ ■ k/(2|q|). The 
crucial difference between ( [Tlj) and ([l3j) is that in the 
latter case the pairing term creates or annihilates pairs 
with total momentum k. A second-order perturbation 
calculation yields the following for the Josephson energy: 



pq 



^p,q+k/27 1 -p, 



-q+k/2A / *A t , 



E p Eq 



9(Eq — Zq 



Er, 



Eq — Zq 



E v 



Eg - Zg) 

V \E a + Z a 



) + e.c.]. (14) 



The expression ( [Ti] ) looks more complicated than Jl2] ) 
because in the FFLO superconductor i) the up and down 
spin quasiparticles are no longer degenerate and ii) the 
ground state is not a vacuum of the quasiparticles. The 
fundamental difference between ( |l4| ) and (|l2|), however, 
lies in the fact that in (^) the two tunneling matrix 
elements forming the product, T p q+k / 2 and T_ Pj _ q+k / 2 , 
are no longer related by time-reversal transformation and 
therefore are no longer the complex conjugate of each 
other. This is due to the fact that the Cooper pairs in 
the FFLO superconductor carry a finite momentum k. 
As a consequence of this, the terms in the summation of 
( |l4| ) have oscillatory phases, and the Josephson effect is 
suppressed. 



Consider geometry a) again. In this case the Josephson 
energy may be recovered by applying a properly chosen 
magnetic field parallel to the planes. In the microscopic 
theory it is convenient to choose the following gauge: 



A = — z(z x B), 



(15) 



A is then parallel to the planes. In this case the gauge 
field does not change the tunneling matrix elements. Let 
us assume the lower layer is at z — 0, while the upper 
layer is at z = d. Thus A/ = so there is no change in 
the BCS superconductor. Also A u = —d(z x B) is a con- 
stant gauge field, which shifts the momenta of individual 
electrons in the upper plane by a constant amount 



5k=^(z x B), 
Tic 



(16) 



i.e., we need to replace all the q's by q -I- 5k, and q's 
by q + <5k, in Eqs. ( |i~3| ) and (|T4|). We emphasize that 
if the upper layer were isolated, such a shift induced by 
a constant A field is a pure gauge transformation that 
has absolutely no physical consequence; the only reason 
that such a parallel magnetic field has physical effects is 
because electrons can tunnel from one layer to the other, 
and thus go through loops that enclose a finite amount 
of magnetic flux. 

In particular, when <5k = — k/2 (as in Eq. @), 
just as in the BCS superconductor, the Cooper pairs in 
the FFLO state carry zero momentum (but still finite 
gauge-invariant center-of-mass velocity*.); in this case the 
Josephson energy takes the form: 



Ej{<t>) 



pq 

0(-E Q 



T T A* A 



E pEq_l 



d(Eq 



- z„ 



E„ + E„ 



- Z„ 



z, 



q-k/2 



E„ 



E„ 



Z„ 



) +C.C.' 



(17) 



Clearly in the above summation we no longer have the 
oscillatory phases, and the Josephson energy and current 
are restored. Again, similar arguments apply for geom- 
etry b) with the conclusions that the FFLO momentum 
must lie along the junction to restore the Josephson ef- 
fect and, if this is the case, the restoration will occur at 
the field given by Eq. (^. 

Discussion. While arriving at the same basic conclu- 
sions, the microscopic and macroscopic derivations do 
not lead to results that are equivalent in all details. This 
is due to the implicit assumption made in the macro- 
scopic derivation that the Josephson coupling is uniform 
throughout the junction. This assumption was not made 
at the microscopic level; if made, it will restrict the form 
of the tunneling matrix elements T P!q . To illustrate this 
consider an infinite junction. In this case translational 



invariance in the junction requires T Py , 



i(p-q)-tj jt 



where tj is an arbitrary translation vector in the junc- 
tion. Implementing this constraint in Eq. (p~2|) gives 
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Ej((f>) oc <5k,o which is also the result of Eq. ([!]) after 
the spatial integration is performed. 

In our discussion we have assumed that the supercon- 
ductors are in the clean limit, so that we can label single 
electron states by their momenta, and the order param- 
eter of the FFLO superconductor has either a single mo- 
mentum or is a linear superposition of components with 
a discrete set of momenta. The presence of disorder will 
change the situation. In this case it is still true that 
the Josephson effect between a BCS superconductor and 
an FFLO superconductor is suppressed, as in the lat- 
ter Cooper pairs are not formed between single electrons 
states that are related by time-reversal transformation. 
However, the recovery of the Josephson effect by apply- 
ing a magnetic field will be decreased, since in this case 
the order parameter in the FFLO state is formed by su- 
perpositions of components with momentum that varies 
continuously; while a magnetic field can only recover the 
contribution to the Josephson energy of a single compo- 
nent with a fixed momentum. Thus experimental studies 
of the effects we predict should be performed in clean 
superconductors. 

Ideally, we would like to apply this theory to FFLO 
superconductors that are stabilized by an internal ex- 
change field, rather than an external magnetic field. This 
is because in the latter case, unless the system is truly 
two-dimensional, the orbital effect of the magnetic field 
that penetrates into the superconductor will induce vor- 
tices, which will complicate the analysis. We note that 
while experimentally the potential candidates for FFLO 
superconductors reported so far are stabilized by an ex- 
ternal field [p|, ^0|p^Jl^ | , recent theoretical studies suggest 
that an FFLO state may form in the ferromagnetic su- 
perconductor RuSr2GdCu208 (T9|| without the benefit of 
an external magnetic field. Despite the disclaimer made 
above, we expect that the above results will be robust in 
the presence of a weak orbital effect, for example in quasi 
2D systems where there is very weak interlayer coupling. 
In this case the orbital effects will be minimal when the 
magnetic field is applied in the layers. 

In conclusion we have shown that while the Joseph- 
son current between a BCS and an FFLO superconduc- 
tor is suppressed, it is restored by applying a magnetic 
field in the junction. The field strength and direction at 
which the Josephson current recovers depend upon the 
momenta of the FFLO order parameter, allowing the mo- 
menta to be measured. The observation of a field restored 
Josephson effect will provide unambiguous evidence for 
an FFLO phase. 
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FIG. 1. Josephson junction geometries considered in this 
article. For geometry a) the magnetic field in the junction is 
applied in the plane of the superconductors; while for geome- 
try b) the magnetic field is applied perpendicular to the plane 
of the superconductors, but confined inside the junction area. 
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